Phase Evolution in Spatial Dark States 
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Adiabatic techniques using multi-level systems have recently been generalised from the optical 
case to settings in atom optics, solid state and even classical electrodynamics. The most well known 
example of these is the so called STIRAP process, which allows transfer of a particle between different 
states with large fidelity. Here we generalise and examine this process for an atomic centre-of-mass 
state with a non-trivial phase distribution and show that even though dark state dynamics can be 
achieved for the atomic density, the phase dynamics will still have to be considered as a dynamical 
process. In particular we show that the combination of adiabatic and non-adiabatic behaviour can 
be used to engineer phase superposition states. 

PACS numbers: 67.85.De, 03.75.Lm, 42.50.Dv 



Studying the wave nature of localised single parti- 
cles allows fundamental questions of quantum mechan- 
ics to be addressed. Recently experimental progress has 
boosted this area and experiments that can control single 
particles have become available in various systems. These 
include neutral atoms in optical lattices [l|-y| or micro- 
scopic dipole traps 0, [B| , electrons in quantum dots [fjj 
and several others. One of the advantages of ultra-cold 
atomic systems is their purity and low-noise environment, 
which makes them well suited for applications in quan- 
tum metrology and information [7j, |8[ . 

Developing a robust toolbox for engineering using the 
laws of quantum mechanics is therefore an important 
challenge. Compared to classical physics, many appli- 
cations require coherent evolution not to be disturbed 
and one common process needed is a mechanism that al- 
lows transfer of particles between different trapping sites. 
Tunneling is such a mechanism, however in its direct ap- 
plication it leads to Rabi oscillations which make con- 
trolled experiments difficult [9|. Recently a new method, 
termed coherent tunneling adiabatic passage (CTAP), 
has been sugg ested [IMl, 

which is analogous to the 
three level techniques of STIRAP in optical physics [3] ■ 
This technique allows for high fidelity transport between 
different trapping sites, with adiabaticity as the only re- 
quirement. 

While STIRAP is a well investigated technique in op- 
tical systems 14|, its translation into the atom optical 



realm offers many new and interesting degrees of free- 
dom to be explored. Recently a number of studies have 
focussed on the effects of non-linear dynamics on the 
transfer process [HI, [HI- Here we add another degree 
of freedom by studying states with non-trivial phase dis- 
tributions and show that the CTAP process is not ro- 
bust with respect to conserving the phase and therefore 
the functional form of the density distribution. However, 
we also show that the process can still be used to con- 
trol the phase of the quantum state in question. Phase 
engineering has become an important technique in the 



area of quantum computing recently and the convenience 
of adiabatic techniques is that if the associated energy 
eigenvalue is zero, one can make use of the usually much 
smaller geometrical phases 1(3- HI • While this is true for 
optical systems, we will show that one has to be more 
careful in atom-optical settings. 

Tunnelling of an individual vortex is an interesting 
problem in a number of systems, including Bose-Einstein 
condensates and Josephson junctions. The escape of a 
single vortex from a pinning potential in a Josephson 
junction has been investigated experimentally [19j and 
recently the tunnelling of a vortex in Bose-Einstein con- 
densate has been the subject of numerical study for dou- 
ble well potentials [20|]. Salgueiro et aJ. found that the 
topological defect is preserved on tunnelling of the con- 
densate and, in fact, can be replicated in such a way that 
each potential minimum has a single vortex (2fjj . 

In the following we will first give a brief introduction 
to the CTAP idea and show that the standard approx- 
imation of a three level system is good for the density 
transport. We will then compare this approximation to 
the full integration of the Schrodinger equation for the 
problem and identify the problems for phase stability. 
Finally we will extend the work to look at systems with 
small non-linearities and conclude. 



COHERENT TRANSPORT 

The CTAP process for cold atoms considers a system 
of three micro-traps, between which a single particle can 
tunnel (see Fig. [IJ . The strength of the tunneling is de- 
termined by the distance and barrier height between the 
individual traps. If one assumes all traps to be of the 
same shape, i.e. have resonant energy levels, and only 
allows for adiabatic processes, one can write the Hamil- 
tonian in terms of the asymptotic eigenstates of the in- 
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FIG. 1: (Color online) Schematic setup for CTAP for an 
atomic state carrying one quantum of angular momentum. 
The atom is initially located in the trap on the left hand side 
and all other traps are considered empty. Tunneling only oc- 
curs between nearest neighbour traps. 
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dividual traps, \L), \M) and \R), as 

Jlm 
H = | Jlm Jmr 
Jmr 



(1) 



The Jij are the tunneling matrix elements between neigh- 
bouring traps and we assume the tunneling between non- 
neighbouring traps to be negligible. The eigenstates and 
eigenvalues of this Hamiltonian are well known (T3 | and 
we will focus here on the so-called dark eigenstate given 
by 



\D) = cosd\L) + sirid\R) 



(2) 



This state only has an indirect contribution from the cen- 
tral trap through the mixing angle, tan# = Jmr/ Jlm, 
and its energy eigenvalue is given by Ed = 0. If one 
allows for time-dependent tunneling rates, this state can 
be used to transfer a particle from, say, trapping site \L) 
to \ R) with large fidelity, even in the presence of noise 
Note that during this transfer the particle has 
no probability to ever be in the state \M). 

To remind the reader, let us briefly review this trans- 
fer process: the tunneling frequencies can become func- 
tions of time either through a time-dependent variation 
of the distance between the individual traps, through a 
modulation of the respective barrier heights or through a 
combination of both of them. Here we will assume that 
the trap positions change in time and, since we will as- 
sume piecewise harmonic potentials, this will also lead to 
a decrease in barrier height. If we therefore first decrease 
the distance cLmr and, with a delay, the distance (Ilm, 
we create the familiar counter-intuitive STIRAP timing 
sequence for the values of the tunneling strengths Jmr. 
and Jlm (see Fig. [5] (a) and (b) and calculations below). 
During this process the mixing angle 8 changes from 
to 7r/2 (see Fig. [2] (c)), which allows a particle initially 
trapped on the left hand side to be transferred to the 
trap on the right hand side. This is the essence of the 
celebrated STIRAP technique. 

To calculate the tunneling frequency between two traps 
as a function of trap distance, let us turn to the exactly 



FIG. 2: (Color online) (a) distance between the left-middle 
and the middle-right trap, (b) tunneling strength between 
neighbouring traps derived from the analytical model de- 
scribed in the text, (c) mixing angle 6 and (d) time-dependent 
eigenvalues of the Hamiltonian (fl) 



solvable model system of piecewise harmonic traps 21 
22j, which also guarantees good approximate resonance 
between the levels at any point in time 
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The eigenfunctions of this potential are known to be given 
by parabolic cylinder functions [23| 



*l> 1 (x) = N 1 D u 
Mx) =N 2 D V 



2mcj , gL 

— {X+ 2\ 
2mw d 

— {x -~2 ] 



for x < , (4) 
for x > , (5) 



where Ni and N 2 are the normalisation constants. The 
eigenenergies are given by E/hui = v + ^ and the 
quantum numbers, are determined by requiring that 
the logarithmic derivatives of the two wavefunctions are 
equal at x = 



x=0 



IJ'2 



(6) 



x=0 



By solving this condition we can exactly calculate the 
tunneling strength between the two traps at any time 
during the adiabatic process. If we then diagonalise 
eq. ([1]) we find the eigenvalues displayed in Fig. [2fd), 
where the dark state with the eigenvalue of zero is clearly 
visible. Since we assume that the whole process is carried 
out adiabatically, the system will be in this state at any 
point in time. 



Angular Momentum State 

In order to write down the Hamiltonian ([T]) for a re- 
alistic system a number of approximations must hold. 
First, the eigenstates of the isolated traps have to be in 
close resonance, and second, the levels of the individual 
traps have to be chosen such that they are well isolated 
from all other available states in the system. The sec- 
ond condition can almost always be fulfiled by making 
the process more adiabatic and the first one translates 
into the simple requirement that all traps have the same 
trapping frequency (in case of harmonic traps). For re- 
alistic traps, however, this is problematic, as potential 
forces often add when they start to overlap. Several so- 
lutions have been proposed for this, inclu ding the use of 
time-dependent compensation potentials [12:] - Here we 
assume that this is experimentally possible, as it allows 
us to isolate the physics relevant to the dynamics of the 
phase. 

To examine the influence of the CTAP dynamics on the 
phase distribution of a quantum state let us first carry 
out a full numerical integration of the Schrodinger equa- 
tion for a general system. This will at the same time 
function as a control mechanism for the approximations 
made above, in particular the fact that we neglected non- 
nearest neighbour coupling and treated tunneling as a 
one-dimensional process. For numerical simplicity, how- 
ever, we will only consider a two-dimensional setup here, 
as this will allow us to capture the main physical pro- 
cesses. In scaled co-ordinates Schrodinger's equation is 
therefore given by 

ij t i>{x, y) = (-^V 2 + ~V(x, y)j ^(x, y) , (7) 

where V{x, y) is the trapping potential of the three traps 
in the linear configuration and which fulfils the conditions 
outlined above. For generality, all lengths are scaled with 
respect to the ground state size of the individual har- 
monic traps, ao = yjh/mui, where m is the mass of the 
particle and uj the trapping frequency of the individual 
harmonic oscillator potentials. All energies are scaled in 
units of the harmonic oscillator energy, Eq = Tioj. 

To be specific, let us choose an initial state that carries 
a single quantum of angular momentum and therefore has 
a phase distribution that increases by 2tt for a closed loop 
around the centre of the state. This state is initially lo- 
cated in the trap on the left hand side. The results of 
the numerical integration are summarised in Fig. [3] and 
show surprising dynamics: while the process still leads 
to a 100% transfer of the probability amplitude to the 
trap on the right hand side (not shown), the angular mo- 
mentum of the final state oscillates continuously between 
clockwise and counter-clockwise depending on the over- 
all duration of the process (see upper part of the figure). 
Four examples of this change in phase and density associ- 
ated with four different durations are shown in the lower 
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FIG. 3: (Color online) Angular momentum as a function of 
the overall duration of the process (upper plot) and the final 
states in the rightmost trap at points A, B, C, and D (lower 
plots). 

half of the figure. In the first example (A) the system is 
in exactly the same state as it started out with, whereas 
in (B) the circulation of the flow has been reversed by the 
CTAP process. The plots (C) and (D) show the situation 
where the final state is in a superposition between clock- 
wise and counter-clockwise rotation, which also leads to 
a strongly modified density distribution. 

This behaviour might seem surprising at first, as angu- 
lar momentum is usually thought of as a conserved quan- 
tity. However, in non-rotationally symmetric geometries 
this conservation law does not hold and in fact can lead 
to interesting dynamics for vortices in anisotropic poten- 
tials (24l - [27j |. As in our example the particle is tunneling 
between rotationally symmetric potentials, it is not a pri- 
ori clear where the necessary asymmetry comes from. 

To gain insight into this behaviour, let us return to 
the three state model and consider the transport of a 
particle in a harmonic oscillator potential carrying one 
unit of angular momentum. In two dimensions its energy 
is given by E = 2Tiuj and the state is doubly degenerate, 
i>n x n v = iho an d ^n, = "001 ■ We therefore have to 
write the most general wavefunction as the superposition 

ip(x,y) = (p 1 (x)ip (y) + i(p Q (x)ip 1 (y)e~ ie , (8) 

where the one-dimensional, single particle eigenfunctions 
of the harmonic oscillator for the ground and first excited 
state are given by ip n . We have also allowed for a relative 
phase 9 between the two degenerate states, however we 
can set this initially to zero without loss of generality. 

As our traps are arranged in a linear configuration 
along the y-axis we can assume that the dynamics in the 
different spatial directions decouple. Tunneling therefore 
needs to be taken into account only for the wavefunction 
part in the y-direction and we find that the Hamiltonian 
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can be split into one for the ground state parts of the 
wavefunction and one for the excited states 
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Unlike in the previous section, we cannot simply set the 
diagonal elements equal to zero as now both the ground 
state and first excited energies are involved. However, 
each of the Hamiltonians has still a dark eigenstate with 
the eigenvalues eo and ei, respectively. If initially a single 
particle is in the trap on the left hand side in the state 
given by eq. (jl]), after the CTAP process, the Hamiltoni- 
ans above lead to 



¥>?(!/) 



vo(v)e- if " f eo{t ' )dt ' 



(11) 
(12) 

(13) 

(14) 



where tf is the overall duration of the process and the 
superscripts L (left) and R (right) refer to the traps the 
wavefunction is localised in. Note that although we have 
used scaled energy units, we reintroduce u into the above 
expressions for clarity. One can immediately see that if eo 
and ei are not independent of time, the wavefunction ac- 
quires a relative phase between the two degenerate states 
that is dependent on the overall duration of the CTAP 
process 

+ i^(x)^(y)e-i^ tf+S ° feiit ' )dt '] . (15) 

For easier understanding let us rewrite the above state 
by defining a global and a relative phase as 



7 = 2 ut f + 



e (t')dt' 



(16) 



e = -wt f + / [e (f) - ei(t')K , (17) 
Jo 

so that we can write the final state as 

4(x, y; t f ) = e*i [^{x)^{y) + i^(x)^(y)e ie ] . 

(18) 

One can clearly see from this that any slight deviation 
in the difference between the asymptotic energy levels of 
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FIG. 4: (Color online) Evolution of eo (blue dashed), ei (green 
dot-dashed) and eo— ei (red solid) over the course of the CTAP 
process. Both energies are initially separated by uj = 1, but 
do not maintain this separation due to the modification of the 
potential as the traps move closer together. 



the individual traps will have a significant effect on the 
final wavefunction. As, in particular, in realistic situa- 
tions the shape of the individual potentials most likely 
strongly changes, one can expect to be unable to control 
the final form of the wavefunction. The CTAP process is 
therefore unstable with respect to states with non-trivial 
phase distributions (28j . However, as our example shows, 
this does not have to be a random process and in fact it 
can be used to engineer the phase of the wavefunction 
deterministically: by slightly changing the overall time 
of the process, one can cycle through all possible states 
for the fixed energy E Q = 2hu>. 

Since in our simulations the potentials are piecewise 
harmonic, it is not immediately clear where the change 
in energies comes from. Let us therefore in the following 
carefully examine our model and determine the influence 
of various other degrees of freedom. The crucial point 
is that during the CTAP process the energy eigenstates 
in the y-direction slightly change, since the potentials we 
are considering are only piecewise harmonic. By bring- 
ing them closer together, the height of the barrier be- 
tween them changes, leading to a different asymptotic 
eigenstate. In fact, the eigenstates are very close to the 
parabolic cylinder functions defined above. As the first 
excited eigenstate is closer to the local maximum sepa- 
rating the traps, it will be shifted differently compared 
to the ground state. While the difference might only be 
small, the integration over a long, adiabatic time interval 
will lead to a large value for the integral. To demonstrate 
this, we show in Fig. 0] the energies for an atom in the 
state ipo(x) and fi{x) during the CTAP process. It can 
be clearly seen that at the time of approach the energy 
eigenvalues slightly changes, as assumed above. In addi- 
tion, we show that their difference, which is the integrand 
of Eq. [T71 is not constant and thus gives rise to the time 
dependant 6(tf). 

From the above argument it follows that the effect 
should not only depend on the overall time of the pro- 
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FIG. 5: (Color online) Final angular momentum as a function 
of time for different minimum distances between the traps 
dmin = 2.0 (full line), d m in = 2.2 (dashed line) and d m i n = 
2.4 (dotted line). The offset of the different curves has been 
shifted for clarity. The inset shows the increase of the period 
with increasing d m i n . 



cess, but also on the minimum distance to which the traps 
approach each other. In Fig.O we show the resulting os- 
cillations of angular momentum for a number of different 
values for d m i n . As the minimum distance is increased, 
the period of the oscillations increases as well (see in- 
set of Fig. [5]), which is in accordance with the explana- 
tion given above: a larger minimum distance between the 
traps means that the energy levels are deviating less from 
the asymptotic values for perfect harmonic potentials. 



Non-linear CTAP 

Let us finally discuss the effect of a non-linearity on 
the evolution of the phase. The paradigm of an atomic 
non-linear system of well defined phase is a Bose-Einstein 
condensate and its dynamics can be described by the so- 
called Gross-Pitaevskii equation 



Hip(x,y) 



V(x,y) + ^\ 2 



ip{x,y). (19) 



Here U is a measure for the non-linearity, which is re- 
lated to the scattering strength between the atoms [29j |. 
While the CTAP process requires resonances of the three 
asymptotic ground state at any time, non-linear samples 
break this due to the time dependence of a chemical po- 
tential, /X, during the tunneling process. For large chem- 
ical potentials the whole process therefore defaults and 
one has to resort to different techniques for restoring the 



resonance during the process 13, 30]. As we are only in- 
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FIG. 6: (Color online) Angular momentum as a function of 
the overall time the CTAP process takes for different non- 
linearities. 



to only small non-linearities (fi <C hui), for which about 
100% transfer can still be reached. 

In Fig. [6] we show the effect a small, but increasing 
non-linearity has on the final phase. The first thing to 
notice is that the oscillatory behaviour of the angular 
momentum does not get immediately suppressed by the 
non-linearity, and that the only effects are an offset and a 
small reduction in periodicity (not shown). We therefore 
speculate that the effect described in this work can also 
be used to create vortex superposition states in Bose- 
Einstein condensates in a controlled way, as long as the 
samples are only weakly interacting. As it is well known 
that vortex oscillations in an anisotropic potential can 
be suppressed for large enough non-linearities 



24, 2fi|, it 



might be interesting to study this effect in the presence 
of resonance restoring techniques. 



Conclusion 

We have shown that the spatial CTAP process for 
atoms in microtraps does not conserve the wavefunction 
form for states with non-trivial phase distribution. This 
is due to an unavoidable small time dependence of the 
asymptotic eigenstates of the individual traps, which re- 
sults from the overlap between them. As any currently 
suggested realistic systems shows large deformations due 
to the overlap, our work directly applies to current ex- 
perimental setups. 

At the same time we have shown that this instability 
behaves deterministically with respect to a change in the 
overall time of the CTAP process and can therefore be 
used to create well defined angular momentum superpo- 
sition states. Such states can have applications in quan- 
tum information and have recently seen a surge in inter- 
~ 32 



est 
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terested in the phase dynamics, we will restrict ourselves 



We have also shown that the presented pro- 
cess even survives in the presence of small non-linearities, 
allowing therefore to superpose vortices in Bose-Einstein 
condensates. 
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